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ABSTRACT 



Aims. We analyse the magnetic field geometry of the magnetic DA white dwarf PG 1015+014 with our Zeeman tomography method. 
^SJ Methods. This study is based on rotation-phase resolved optical flux and circular polarization spectra of PG 1015+014 obtained with FORSl 
at the ESO VUT. Our tomographic code makes use of an extensive database of pre-computed Zeeman spectra. The general approach has been 
described in Papers I and II of this series. 

Results. The surface field strength distributions for all rotational phases of PG 1015+014 are characterised by a strong peak at 70 MG. A 
separate peak at 80 MG is seen for about one third of the rotation cycle. Significant contributions to the Zeeman features arise from regions 
O ^ with field strengths between 50 and 90 MG. We obtain equally good simultaneous fits to the observations, collected in five phase bins, for 
two different field parametrizations: (i) a superposition of individually tilted and off-centred zonal multipole components; and (ii) a truncated 
multipole expansion up to degree 1 = 4 including all zonal and tesseral components. The magnetic fields generated by both parametrizations 
exhibit a similar global structure of the absolute surface field values, but differ considerably in the topology of the field lines. An effective 
photospheric temperature of T^s = 10 000 ± 1000 K was found. 

Conclusions. Remaining discrepancies between the observations and our best-fit models suggest that additional small-scale structure of the 
magnetic field exists which our field models are unable to cover due to the restricted number of free parameters. 
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1 . Introduction challenged by the recent detections of kilogauss-size fields in 

several MWDs as well as in their direc t progeny (central stars 
In about 170 of the 5448 white dwarfs (WDs) listed in pf planetary nebula e and h ot subdwarfs.lAznar Cuadrado et all 
the Web Version' of the ViUanova White Dwarf Catalog EoO?, -Jordan et all Eool lO'Toole et al.L l2005l) . Undoubtedly, 
magnetic fields between 2kG-1000MG have been de- further theoretical and observational efforts are required in or- 
tected, corresponding to a fraction of - 3 % ( McCook & Sion ^j^^^ ^^^^ j^gj^^ ^„ ^1^^ ^.^j^ j^agnetic fields play in the 

1999;.Wickramasinghe & Ferrariq,200Q^.Vanlandingham et al, ^^^^^^ post-main sequence evolution. For the present 



20051). Most of the WDs have not been scrutini singly exam- purpose, we consider MWDs as stars displaying a field strength 

ined for the presence of a magnetic field, however, and a sta- g > i y[Q 

tistical study suggests that the true fractional incidence could „ , . . . j, . r^^rr^ n ^ ^ 

... ^^Z. ITT-; nrm^ rrr — ~, nr;„„^ Due to the intrinsic taintness of WDs, 8-m class telescopes 

be as high as 20% jLiebert et alJl200a ISchmidt et alJ l2003). ... . , .... ,.' , , • I- 

^, ... . ' — , ,„,Tr~ — . , , — ; — ~ — ^ , , are required in order to record high-quality spectropolarimetric 

The magnetic white dwarfs (MWDs) are widely believed to be , ^ rr- ■ , , ^- t. ■ c ^ j- 

. . data with sufhcient time resolution as a basis for studies of the 



the successors of the chemically peculiar magnetic Ap stars, 
which are the only main sequence stars to show substantial 
globally organised magnetic fields. However, this scenario is 



magnetic field geometry. In the course of our Zeeman tomogra- 
phy programme we have conducted observations for a number 
of isolated (non-accreting) and accreting MWDs at the ESO 

"~ ; ~ ^ ^ , ., VLT with FORSl in the spectropolarimetric mode. 

Send offprint requests to: F. Euchner, e-mail: 

feuchnertostro. physik.uni-goettingen.de I" the present work - the third paper of our series on 

* Based on observations collected at the European Southern Zeeman tomography - we present an application of our code 

Observatory, Paranal, Chile, under programme ID 63.P-0003(A). to observations of the non-accreting white dwarf PG 1015+014 

' http://www.astronomy.villanova.edu/WDCatalog/index.html, and find a field geometry that deviates strongly from cen- 

January 2006. tred dipoles or quadrupoles. In the first paper of the series, 
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Fig. 1. Left panel: Flux spectra of PG 1015+014 from May 1999. The uppermost 1 1 spectra have been shifted upwards by 2 flux 
units, respectively. Right panel: Combined flux and circular polarization spectra of PG 1015+014 from the May 1999 observa- 
tions. These spectra, which have been collected into five phase bins, will be used as input spectra for the Zeeman tomography 
procedure. For clarity, the uppermost four flux spectra have been shifted upwards by 1.5 flux units, respectively. The tick marks 
below the top axis indicate the wavelengths that were used to adjust the continuum flux level of the the model spectra to the 
observations. 



we have demonstrated that Zeeman tomography is a suitable 
method to recove r field geometries by analysing synthetically 
generated spectra jEuchner et alJ20(M hereafter Paper I). In a 
first application of this theory, we have derived a quadrupole- 
dominated field st ructure with a preva iling field of ^ 16 MG for 
HE 1045-0908 ( Euchner et alJ2005l hereafter Paper II). 

The magnetic DA white dwarf PG 1015+0 14, discov- 
ered in the Palomar Green survev (iGreen et alJ [T987). was 
observed by IWickramasinghe & CroppeJ ( 1 19881 hereafter 
WC88) with the RGO spectrograph at the AAT in the wave- 
length range 4000-7000 A. Their phase-resolved spectroscopy 
and low-resolution circular polarimetry revealed significant 
modulations in flux and circular polarization (V//) over the 
rotation cycle. From nearly sinusoidal oscillations of the 
wavelength-averaged degree of circular polarization between 
-1.5% and 1.5% the authors derived a rotational period of 
f rot = 9 8.7 min, which was later confirmed with higher accu- 
racy by Sch midt & Norsworthvl ( Il99ll) who used white-light 
circular polarimetry. In the individual polarization spectra of 
WC88, is ^ 5 % in the continuum and up to ^ 10 % in in- 
dividual features. They fitted theoretical MWD model spectra 
to the observations and found an obliquely rotating magnetic 
dipole model with a polar field strength of B^^^ = 120 + 10 MG 



and an almost equator-on view to be the best-fitting field geom- 
etry. Remaining discrepancies between observations and model 
spectra were attributed to higher-field regions superimposed on 
the dipolar field structure. Our analysis provides a substantially 
improved insight into the field structure of PG 1015+014. 

2. Observations 

We have obtained spin-phase resolved circular spectropo- 
larimetry of PG 1015+014 with FORSl at the ESQ VLT on 
May 15, 1999. The spectrograph was operated in spectropolari- 
metric (PMOS) mode, with the GRIS_300V+10 grism and an 
order separation filter GG 375, yielding a usable wavelength 
range -3850-7900 A. With a slit width of 1" the FWHM 
spectral resolution was 13 A at 5500 A. The observational data 
have been reduced according to standard procedures (bias, 
flat field, night sky subtraction, wavelength calibration, atmo- 
spheric extinction, flux calibration) using the context MOS of 
the ESO MIDAS package. The instrumental setup and the data 
reduction are analogous to those employed for our analysis of 
HE 1045-0908 (Paper II). 

We secured a sequence of 12 exposures with an exposure 
time of 8 min each, covering a full rotation cycle. We were able 
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Fig. 2. Zeeman maximum (0 = 0.25, te/f panel) and minimum (0 = 0.86, right panel) flux (/j) and circular polarization (V//) 
spectra for PG 1015+014, plotted along with the theoretically predicted field-dependent transition wavelengths for Balmer ab- 
sorption lines {A-B curves). Transitions corresponding formally to |A/| = 0, 2, 4, . . . with the zero-field angular momentum quan- 
tum number / have been plotted as dashed curves in the bottom parts of the figures. Filled circles denote unambiguous identifica- 
tions of transitions, open circles less certain ones. In the top parts, theoretical circular polarization spectra are shown for compar- 
ison (shifted upwards for clarity, see the text for the adopted model parameters), is given in units of 10"'^ erg cm"^ s"' A"' . 



to reach a signal-to-noise ratio S IN ^ 100 for the individual 
flux spectra (Fig.^ left panel). The wavelength-dependent de- 
gree of circular polarization ( V/ /) was computed from two con- 
secutive exposures - with the retarder plate positions differing 
by 90° - in order to eliminate Stokes parameter crosstalk, thus 
yielding six independent phase bins. The flux spectra do not 
show the typical absorption signature of low- and intermediate- 
field MWDs with the Ha n component at the zero-field wave- 
length surrounded by broader cr" and cr+ troughs, but exhibit 
a variety of distinct sharp lines scattered over the whole op- 
tical range which vary in strength, position, and shape with 
the rotational phase. The continuum circular polarization dif- 
fers significantly from zero for most phases. Both phenomena 
are characteristic of a high-field o bject with B t 50 MG (e.g. 
IWickramasinghe & Ferrariol200(]t) . 



We collect the observed flux and circular polarization spec- 
tra into five nearly-equidistant phase bins (0 - 0.05, 0.25, 
0.45, 0.66, 0.86). The phas es refer to the ephemeris of 
ISchmidt & Norsworflivl (1 1 99 ih with (p -Q denoting the change 



from positive to negative continuum polarization^. These flux 
and polarization spectra (Fig. ^ right panel) form the ba- 
sis of our tomographic analysis. We investigate two particular 
phases in greater detail: = 0.25, when the Zeeman features 
are strongest ("Zeeman maximum"); and <p - 0.86, when they 
are weakest ("Zeeman minimum"). 

3. Qualitative analysis of the magnetic field 
geometry 

As a first step of the analysis, we compare the positions and 
strengths of the most prominent features in the observed spec- 
tra with the expected field-dependent wavelengths of the hy- 
drogen tra nsitions A (B), henceforth referred to as A-B curves 
(fForster et alJl984 tRosner et al. 1984; Wunner et al. 1985). In 
Fig. 121 we show the observed flux and circular polarization 
spectra along with the A-B curves for (p - 0.25 (left panel) and 
0.86 (right panel). Transitions that could be unambiguously 
identified with specific spectral features have been marked with 

- These authors used a different convention for the sign of the cir- 
cular polarization than we did, since in their data V// changes from 
negative to positive values at = 0. 
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Table 1. Positions of observed Zeeman features and corre- 
sponding Balmer transitions. A colon denotes a less certain 
identification. 



A I K Transition 



A I K Transition 



(0 = 0, 
4610 

4805 
5117 
5326 

5475 
6193 
6268 
6326 
6927 
7085 
7369 
7594 



25j 
H/3(2so-; 
He(2p+,- 
H/?(2po- 
H/?(2so-H 
Ha (2p_i- 
H0(2p+i- 
Ha (2po- 
Htt (2p±i - 
Ha (2s 0- 
Ha (2po- 
Ha (2p+i ■ 
Ha (2s 0- 
Ha (2p_i- 
Ha (2po- 



4fo) 
■.4so) 
4d_i) 
4f-i) 
-.3d„) 
-^4do) 
•3so): 
-.3d±,) 

3po) 
'3do) 
-»3so) 
3p-i) 
-» 3d_2) 
•3d_i) 



(0 = 0.86j 

4755 H8(2po^4d_i) 

5115 H/3(2so-»4f_i) 

5480 HQ'(2p()^3so): 

6088 HQ'(2p±i-»3d±i) 

6198 Ha(2so-»3po) 

6243 Ha(2po^3do) 

6864 HQ-(2p+i -»3so) 

7060 Hq'(2so^3p_i) 

7405 HQ'(2p_i^3d_2) 

7705 HQ'(2po-^3d_i) 



filled grey circles and are listed in Table[T] Possible or unlikely 
identifications are displayed as open circles. In the presence 
of an electric field additional transitions can occur which for- 
mally correspond to |A/| = 0,2,4, . . ., with I being the angular 
momentum quantum number in the zero-field nomenclature. 
These components have been plotted as dashed curves in the 
bottom parts of Fig.|2l 

In the Zeeman maximum spectrum (0 = 0.25), three Ha n 
and four He transitions with A > 6000 A are clearly visi- 
ble as positive peaks superimposed on a negative polarization 
continuum. In particular, the narrow feature at 7085 A which 
corresponds to a local extremum of the 2s()^ 3p_i A-B curve 
allows a reliable estimation of the prevailing field strength. For 
A < 6000 A, we identify three strong features showing positive 
polarization (at 4610, 4805, and 5117 A) with H/3o-+ transi- 
tions. Our identifications of several line features differ from 
those suggested by WC88 for the same rotational phase. The 
feature at 5326 A shows a noticeable negative polarization and 
is interpreted as a blend of a Ha cr" and a Hp cr+ component, 
with the 2p_i^3d() transition dominating. The strong line 
at 5475 A with positive polarization probably arises from the 
Ha n transition 2po^ 3so. The feature at 5930 A reported by 
WC88 is not present in our data. Furthermore, there is a distinct 
feature with a positive peak in polarization at 5200 A which 
they attribute to the Ha cr" transition 2p+i 3d+2- We ques- 
tion this identification because the Ha cr" components pro- 
duce lines with strong negative polarization. However, no ob- 
vious match with a specific A-B curve can be found in Fig.|2l 
so the origin of this feature remains unclear We do not at- 
tempt to identify spectral features at wavelengths A < 4600 A 
because the number of candidate transitions from the overlap- 
ping Hj6 and Hy manifolds is too large. Taking into account all 
these identifications, we conclude that the distribution of the 
corresponding field strengths must be centred quite sharply at 
^70MG. To illustrate this, we show in the top part of Fig.|2 
(left panel) a theoretical circular polarization spectrum for a 



Table 3. Best-fit magnetic parameters for a superposition of 
three individually off-centred, non-aligned dipoles (labelled 
DDD sum in Fig.O. The individual dipole components are de- 
noted by Dl, D2, and D3. The uncertainties in the last digit 
are given by the values in brackets. The best-fit inclination is 
23°±4°. 





Dl 


D2 


D3 


Ko> (MG) 


-40(2) 


92(5) 


-38(3) 


0^ n 


44(4) 


63 (2) 


63 (5) 


n 


339(2) 


276 (6) 


134(3) 


x'oS (Rwd) 


0.04(1) 


-0.012(2) 


0.27(1) 


y'os (^wd) 


0.35(3) 


-0.136(8) 


0.080(7) 


Z off (^Wd) 


0.33(1) 


-0.28 (3) 


0.21 (2) 



single value B = 69MG with lA = 51°, and T^s = 10 000 K, 
where i// denotes the viewing angle between the magnetic field 
direction and the line of sight. 

In the Zeeman minimum spectrum (0 - 0.86), the absorp- 
tion signature for A > 6000 A is similar to that at Zeeman max- 
imum, although the Ha n line components are not visible as 
separate peaks in V/I, but rather produce a broad depression of 
the overall positive continuum polarization. The Ha cr+ com- 
ponents are wider than those at cp = 0.25. Both effects suggest 
a broader distribution of field strengths at Zeeman minimum. 
A theoretical circular polarization spectrum for a single value 
B = 82 MG, (A = 129°, and T^ti = 10 000 K provides a good in- 
dication for the prevailing field strength and direction. As can 
be seen from the change of polarity in the continuum circular 
polarization, the net magnetic field direction has changed sign 
along the line of sight. The fact that for both phases, which 
are separated by about half a rotation cycle (A0 = 0.39), the 
field distributions seem to be clearly peaked at values differing 
by only ^ 20 % indicates a field geometry more complex than 
centred or off-centred dipoles or quadrupoles. 

4. Zeeman tomography of the magnetic field 

Our Zeeman tomography synthesises the observed spectra in a 
best-fit approach. It makes use of an extensive pre-computed 
database of theoretical flux and circular polarization spectra 
of magnetic white dwarf atmospheres, with B, i/r, Teff, \ogg, 
and the direction cosine fi - cos § as free parameters, where 
§ denotes the angle between the normal to the surface and the 
line of sight. The three-dimensional grid of 46 800 Stokes / 
and V profiles covers 400 B values ( 1-400 MG, in IMG 
steps), nine ifr values (equidistant in cos i/r), and 13 temperatures 
(8000-50 000 K) for fixed logg = 8 and = 1. Limb darken- 
ing is accounted for in an approximate way by the linear inter- 
polation law Ia{m)/^a,ii=i - a + b^ with constant coefficients 
a - 0.53 and b - 0.47. The best fit to the absolute flux distri- 
bution of PG 1015+014 was obtained with an effective tem- 
perature of Teff - 10000 + lOOOK. We adopt this value in the 
subsequent analysis. 

In order to determine the misfit between the observations 
and the model spectra for a given set of magnetic parame- 
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Table 2. Best-fit magnetic parameters for the different parametrizations of the magnetic field shown in Fig. 13 The uncertainties 
in the last digit are denoted by the values in brackets. For each model, a short description is given in the first line. The best-fit 
inclinations are 24° ± 3 ° (model 1), 45° + 3 ° (model 2), and 36° + 4 ° (model 3). The superscripts 'd', 'q', and 'o' refer to dipole, 
quadrupole, and octupole, the subscript 'pol' to the polar field strength. The last three columns give the offsets in units of the 
white dwarf radius. Formulae for the field structure are given in Paper I. 





0d 






01 






0° 


o° 


^ off 


/off 


Z'off 


(MG) 


n 


(°) 


(MG) 


n 


o 


(MG) 


(°) 


(°) 


(^wd) 


(^wd) 


(^wd) 



(1) D offs ( off-centred dipole) 

-97(2) 85(3) 77(4) - - - - - - -0.113(5) -0.0036(2) 0.162(4) 



(2) DQO na ctrd (non-aligned, centred combination of dipole, quadrupole, and octupole) 

-1.4(1) 53(13) 325(6) 13.9(6) 26(2) 115(5) 174(1) 65(3) 94(5) _ _ _ 

(3) DQO na offs (non-aligned, off-centred combination of dipole, quadrupole, and octupole) 

-38(2) 85(3) 82(5) -15(3) 85(5) 21 (3) 171 (6) 75(6) 102(5) 0.060(5) 0.011(1) 0.081 (4) 



Table 4. Best-fit magnetic parameters for a truncated multipole 
expansion up to degree / - 4 (labelled M in Fig.|3}. The coef- 
ficients and are in MG. The best-fit inclination is 47°. 
The tilt and the azimuth of the multipole axis relative to the 
rotational axis are 22° and 191°. 
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0.6 
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ters, we use the classical reduced as our penalty function. 
The optimisation problem of finding the best-fitting parameters 
is solved using an evolutionary strategy, implemented by the 
evoC library (Trint & Utecht 1994)^. We assign equal statistical 
weight to the flux and polarization spectra for all five rotational 
phases and to all wavelengths within the individual spectra. The 
statistical noise of the observations entering the;^'^ function has 
been estimated by comparing the observed spectra after the 
appUcation of a Savitzky-Golay filter of 20 A width with the 
unprocessed versions. The flux level of the model spectra has 
been adjusted to the observations at the wavelengths marked 
with ticks at the top of ¥\g.^(right panel). This is an attempt 
to remove differences between observed and model flux spec- 
tra with wavelengths 100 A, as expected from the remain- 
ing uncertainties in the detector response function. We do not 
correct either for differences between the spectra on shorter 
wavelengths, nor do we apply any correction to the polariza- 
tion spectra. The resulting x^^^ values are of the order of 80 
indicating a gross underestimation of the errors which enter 
the computation. Larger observational errors and remain- 
ing systematic uncertainties in the model spectra probably both 

^ ftp://biobio.bionik.tu-berlin.de/pub/software/evoC/ 



contribute to this result. We have experimented with different 
statistical weights assigned to different wavelength regions - 
as specific line features or subsections of the continuum - but 
found no convincing way to better define the goodness of fit. 
We compromise on using the formal x^^^, even if large, as a 
guide line and decide by an admittedly subjective 'by eye' pro- 
cess which individual of similarly good fits to accept. 

4.1. Field parametrization 

We have embarked on two different strategies in order to find an 
adequate parametrization of the magnetic field geometry (see 
also Paper I). The classical approach is that of an expansion 
of the scalar magnetic potential into a series of spherical har- 
monics, characterised by degree I and order m with in = to / 
for each I. Given I and m, two free parameters and /i™ have 
to be specified, o r one parame ter only for the zonal compo- 
nents with m - llLangel Pl987). The approach is powerful, but 
is limited to the description of rather simple structures if one 
truncates the expansion at low values of I. Such a truncation 
is necessary in order to avoid convergence problems of the fit, 
given the rapid growth in the number of fit parameters which 
increases as 1(1 -H 2) for the full I, m-expansion. Conceptually 
simple structures, such as the sum of a quadrupole and octupole 
with their axes inclined relative to each other, cannot be re- 
alised if the series is truncated at low / > 3. As an alternative 
approach, we adopt a hybrid model consisting only of zonal 
(in = 0) harmonics with independent tilt angles and off-centre 
shifts. All tesseral components (m + 0) are ignored in this case. 
Examples are, e.g., a combination of dipole, quadrupole and 
octupole, and also a combination of three dipoles indepen- 
dently inclined and offset from the centre. Further details on 
the field parametrizations are given in Papers I and II. 

In the following Section, we proceed systematically from 
simple structures, as centred or off-centred single zonal com- 
ponents, over the sum of such components, individually tilted 
and/or offset, to a full multipole expansion truncated at Z = 4. 
The number of free parameters varies between 4 and 27. In the 
case of the hybrid model, the fit parameters for each component 
are the polar field strength Bpoi, the tilt © and azimuth <1) of the 
magnetic axis, and, if appUcable, the offsets x'off, y' os, and z'oir 
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Fig. 4. Observed spectra of PG 1015+014 (grey curves) and best-fit synthetic spectra (black curves) using a superposition of three 
individually off-centred, non-aligned dipoles. The upper four flux (circular polarization) spectra have been shifted for clarity by 
2, 4, 6, and 8 (0.1, 0.2, 0.3, and 0.4) units in /i (V/I). 



from the centre plus the inclination of the rotational axis rela- 
tive to the line of sight. For the truncated multipole model, we 
fit the 1(1 + 2) parameters of the model, the two angles of the 
reference axis, and the inclination. 



4.2. Results 

For single centred dipoles, quadrupoles, or octupoles no 
satisfactory fits could be obtained. Nevertheless, we note 
that the best fit with a centred dipole to the single phase 



<p = 0.25 yielded a value of BjJ^, = 131 + 1 MG and an an- 
gle of a = 83° between the magnetic dipole axis and the 
Une of sight, which is co mpatible with the values quoted by 
IWickramasinghe & Cropper (.1988 ). The model is far from op- 
timal, however, in particular for the other phases. 

In Fig.|3lwe show the observed and best-fitting model cir- 
cular polarizations at = 0.25 and 0.86 for selected hybrid and 
multipole parametrizations with increasing complexity. The 
simplest configuration is that of an offset dipole (denoted by 
D offs in Fig.E). The corresponding values of the magnetic pa- 
rameters are listed in Table |2] The overall shape of the con- 
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Fig. 5. Zeeman tomographic analysis of the magnetic field configuration of PG 1015+014 using a superposition of three individ- 
ually off-centred, non-aligned dipoles. Left panel: B-tfr diagram, right panel: absolute value of the surface magnetic field, cosine 
of the angle ij/ between the magnetic field direction and the line of sight, and maximum radial distance reached by field lines in 
units of the white dwarf radius. 



tinuum polarization is reproduced fairly well for all phases, 
but the model fails seriously in the spectral lines, especially 
for (f) - 0.86. For an offset quadrupole, the result was compara- 
bly poor. An offset octupole fared slightly better, but still was 
considered unsatisfactory. While acceptable fits could be ob- 
tained with offset dipoles and quadrupoles for HE 1045-0908 
(Paper II), the magnetic field geometry of PG 1015+014 seems 
to be significantly more complex. 



As a next step, we proceed to the same hybrid field mod- 
els that yielded successful fits for HE 1045-0908. These are 
produced by the superposition of dipole, quadrupole, and oc- 
tupole components and the introduction of a common off-centre 
shift. For the nonaligned centred dipole-quadrupole-octupole 
combination (DQO na ctrd in Fig. |3}, the polarization spec- 
trum at <p - 0.86 is reproduced remarkably well, in particular 
the broad negative dips in polarization at 43 10 and 4740 A, and 
the 2so — > 3p_i transition at 7085 A. The frequency distribution 
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Fig. 6. Observed spectra of PG 1015+014 (grey curves) and best-fit synthetic spectra (black curves) for a truncated multipole 
expansion up to degree I = 4. For clarity, the upper four flux (polarization) spectra have been shifted upwards by 2, 4, 6, and 8 
(0.1, 0.2, 0.3, and 0.4) units in /, (V/I). 



of magnetic field strengths and directions (B-i// diagram, not 
shown) is double-peaked at fields of 69 and 8 1 MG in accor- 
dance with the considerations of Sect.|3l The essence, however, 
is that the reproduction is poor for the other phases, including 
(p - 0.25, although (p - 0.86 fits well. Introducing a common 
offset to the model (DQO na offs in Fig.O improves the fit at 
(f> = 0.25, but shows significant deviations at (p - 0.86. As can 
be seen from the best-fit magnetic parameters in lines (2) and 
(3) of Table |2j both the centred and the shifted hybrid model 
are dominated by the octupole component, again suggesting a 
rather complex geometry. We conclude that the models consid- 



ered so far still do not suffice to achieve a fair reproduction at 
all phases. 

In an attempt to improve the fits, we stayed with hybrid 
models and tried to model a star with several field concentra- 
tions by adding three dipoles which are individually tilted and 
off-centred. The model fares surprisingly well despite the ab- 
sence of quadrupole and octupole components. We conclude 
that the ability to place and adjust three spots meets reality 
rather closely. The fit results for this 'triple dipole' model are 
displayed in Fig.|3 Fig.|5] and Tableland are discussed below. 
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Fig. 7. Zeeman tomographic analysis of the magnetic field configuration of PG 1015+014 using a truncated multipole expansion 
up to degree I = 4. Left panel: B-i// diagram, right panel: absolute value of the surface magnetic field, cosine of the angle if/ 
between the magnetic field direction and the line of sight, and maximum radial distance reached by field lines in units of the 
white dwarf radius. 



Finally, we considered a truncated multipole expansion up 
to Z = 4 which includes all orders with indices m -Q to I. 
This 'truncated multipole' model achieves a y^^^ value simi- 
larly good as the triple dipole model and both provide substan- 
tially better fits than the parametrizations discussed above if all 
rotational phases are considered. In particular, the wavelength- 
dependent level of the continuum circular polarization is repro- 
duced more accurately. The results for the truncated multipole 
model are shown in Figs.|6land0and in Tabled] 



Although both magnetic field topologies are generated by 
entirely different parametrizations (Tables |3] and 0}, the over- 
all appearance of their B-ili diagrams is similar For both con- 
figurations, the visible part of the surface field is dominated 
by an extended area with a rather small variation of the field 
strength that leads to the pronounced peaks at ^ 70 MG in the 
B-if/ diagrams. The phases (p - 0.25 and 0.45 are almost en- 
tirely dominated by this region of constant (although not ho- 
mogenous) field. For the remaining phases, the low- and high- 
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Fig. 8. Best-fit results for individual single phase fits to the Zeeman maximum (0 = 0.25, bottom panels) and minimum (0 = 0.86, 
top panels) phases. The field has been parametrized by a truncated multipole expansion up to degree 1-5 for - 0.25, and up 
to / = 4 for (p - 0.86. Left panels: B-ijj diagram, right panels: absolute value of the surface magnetic field, cosine of the angle 
if/ between the magnetic field direction and the line of sight, and maximum radial distance reached by field lines in units of the 
white dwarf radius. 



field tails of the field strength distribution become more promi- 
nent. In the triple dipole model, the regions with low and high 
fields become manifest as small spots, one of which reaches up 
to ^ 160 MG. A second high-field region appears at the stellar 
limb for phases = 0.66 and 0.86. For the truncated multipole 
model, a similarly prominent high-field spot is not seen, but a 
low-field spot appears close to the stellar limb for (p - 0.66 and 
0.86 and finds its counterpart in the hybrid model. The high- 
field regions are divided into two small spots with field values 
up to ^ 90 MG and pronounced negative cos i/r, leading to a dis- 
tinctive signature in the B-if/ diagram, and two narrow areas at 
the stellar Umb with B up to ^ 150 MG and 245 MG, which be- 
long to high-field spots on the hidden part of the stellar surface. 
This comparison of the triple dipole model and the truncated 
multipole model emphasises the fact that the principal infor- 
mation on the field structure is contained in the B-if/ diagrams. 
These diagrams indicate that the distributions of the magnetic 
field are similar, a conclusion which is impossible to draw from 
the numerical descriptions in Tables|3and|3 

For both configurations substantial deviations from the ob- 
servational data remain, e.g., at the broad dips in V// at 4310 
and 4740 A for phases 4> - 0.66 and 0.86 which both models 
cannot reproduce correctly. While it is quite likely that a larger 
number of free parameters - and, hence, a still further increased 
complexity of the field - will improve the fits, our present op- 
timisation procedure cannot handle more free parameters and 
prevents us to pursue this possibility further. As a consequence. 



it remains unclear what part of the remaining discrepancies, if 
any, may be due to systematic errors in the model spectra. 

A possible way to answer the last question at least in part is 
to perform fits to the flux and polarization spectra for a single 
phase and abandon the requirement that the model should si- 
multaneously fit the other phases. The resulting model may be 
wrong in a global sense, but will provide a more accurate de- 
scription of the field distribution over the visible hemisphere at 
the selected phase. Figure|Slshows the resulting B-ijj diagrams 
at = 0.25 and = 0.86 using the truncated multipole model 
with I - 5 and I - 4, respectively. The corresponding po- 
larization spectra are shown at the bottom of Fig. |3] The 
B-if/ diagrams differ, in fact, from the distributions for the same 
phases obtained from the simultaneous fits to all five phases in 
Figs.|5land0 The most obvious new feature is the appearance 
of a second field maximum at 82 MG for (p - 0.86. At = 0.25, 
the distributions of the field strengths look similar, but the i^- 
distributions and, hence, the field geometries differ (see Figs.|5] 
and|S}. We conclude that a major fraction of the misfits still 
present in our two best global fits is due to the disability of the 
models to appropriately account for the complexi ty of the field . 
The parameter-free spectral synthesis of Dona ti et alJ 
would allow to optimise the fit to a single phase still further 
and answer the quest for the best possible fit with the present 
database, at the expense of a global field solution, however (see 
also Section 5). 

We conclude that the field models used by us are barely 
sufficient to describe a single phase of the observations of 
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Fig. 3. Observed circular polarization spectra for phases 
(f> - 0.25 and 0.86 (grey curves) and best-fit synthetic spectra 
(black curves) for different parametrizations of the magnetic 
field geometry. From top to bottom: off-centred dipole (D offs); 
centred, non-aligned combination of dipole, quadrupole, and 
octupole (DQO na ctrd); off-centred, non-aligned combina- 
tion of dipole, quadrupole, and octupole (DQO na offs); super- 
position of three individually off-centred, non-aligned dipoles 
(DDD sum); truncated multipole expansion up to degree / - 4. 
The lowermost two curves show individual fits to single phases 
with a truncated multipole expansion up to Z - 5 {(p - 0.25) 
and I — 4 (0 = 0.86). All curves except for the bottom one 
have been shifted vertically by suitable amounts in V//, with 
the horizontal dashed lines indicating the respective levels of 
zero polarization. 

PG 1015-1-014, and are certainly not complex enough to de- 
scribe the global field configuration by fits to all five phases. 

Several potential sources of deviations between obser- 
vations and synthetic model spectra have been proposed in 
Paper II and arise both from remaining uncertainties on the 
observational (flat-fielding, flux calibration) and on the theo- 
retical side. Since PG 1015-1-014 is dominated by substantially 
higher magnetic field strengths than HE 1045-0908, additional 
error sources for the theoretical model spectra have to be con- 
sidered. With growing magnetic field strength, for instance, the 
influence of electric fields on line positions and strengths be- 
comes increasingly important. For the case of arbitrarily ori- 
ented electric and magnetic fields no discrete symmetry is left, 
leading to slightly different transition wavelengths and oscil- 



lator strengths than those computed for the diamagnetic case, 
and even to the occurrence of additional dipole transitions 
jFaBbinder & Schweizei|[l995 iBurleieh et alJll99 9<). Lines of 
such type have not been included in the atomic data tables used 
for the computation of our synthetic model spectra and, con- 
sequently, cannot be reproduced by our fits. They are possibly 
responsible for the sharp line features at 5200 A and 5475 A 
and the washed out feature at 5750 A, which are not explained 
by our models (Figs.|2and|6j. 

Another theoretical uncertainty arises from the simpli- 
fied treatment of the field-dep endent bound- free and free- 
free transitions as described in Jordan ('1992"). For the case 
of Grw -i-70°8247 with B ^ 320 MG, Jordan & Merani ( 1995) 
have shown that a more consistent but numericaUy extremely 
expensive treatment of the continuum opacities can yield 
slightly different results for the polarization. We regard the 
uncertainties arising from both potential error sources as rel- 
atively small at fields of ^ 80 MG, but slight deviations from 
our synthetic spectra cannot be ruled out. 



5. Discussion 

In this work, we have analysed high-quality spectropolarimet- 
ric data of PG 1015-(-014 covering a whole rotational period 
with the Zeeman tomography code described in Papers I and 
II. We have achieved good fits to the observations, but require 
a magnetic field geometry that is significantly more complex 
than the popular assumption of centred or moderately offset 
dipoles and quadrupoles propo sed by several authors in the past 
jWickramasinghe & Ferrariab oOO. and references therein). In 
fact, the magnetic field structure of PG 1015-H014 is even more 
complex than that derived for HE 1045-0908 (Paper II) and 
can be successfully modelled by two different parametriza- 
tions: On the one hand, we used a hybrid model of off-centred, 
tilted zonal (m - 0) components, and on the other hand, a trun- 
cated multipole expansion including all mi^Q components. 
Hence, our results further confirm the evidence that the mag- 
netic field structures of MWDs are non-trivial and require 
higher multipole components for an accurate description. We 
have shown, in particular, that a simple oblique dipole model 
as devised by^Wickramasinghe & Cropper ( 1988) does not suf- 
fice to describe the complex Zeeman absorption features. 

Our findings for PG 1015-1-014 show the difficulties that are 
inherent to the description of a star's complex magnetic field 
geometry with only a few numeric parameters. For this ob- 
ject, an appropriate description requires to quote at least the 
range of field strengths that contributes effectively to the spec- 
tral shape, 50 to 90 MG, and the phase-dependent maxima of 
the field strength distributions at 69 and 82 MG. 

For both parametrization strategies the number of free pa- 
rameters required for a good fit approaches the maximum num- 
ber our tomography algorithm can currently handle, and we 
could reach equivalent quality levels for the best fits. The fits 
are sufficiently good to be certain that the best-fitting field ge- 
ometry comes reasonably close to reality, but we could not 
reach the same high quality of the fits as for HE 1045-0908 
(Paper II). 
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We propose two explanations for the remaining differences 
between the observations and our best-fitting models: (i) sys- 
tematic uncertainties in the model spectra arising already in the 
field regime at B ^ 50 MG, as discussed at the end of Sect. 14. 21 
(ii) insufficient spatial resolution of the magnetic field distri- 
bution provided by our field geometry models. Given the lim- 
ited number of free field parameters and the corresponding 
limitations regarding the attainable level of complexity in the 
B-ili diagrams, residuals caused by both effects cannot be dis- 
entangled. However, the fact that fits for a single phase were 
clearly superior to the simultaneous fits for all five phases sug- 
gests that there is room for improvement on the side of finer- 
grained field distributions. Hence, it would be necessary to ex- 
amine the quality of fits to all phases with an increased number 
of free parameters before a reliable estimation of the effects of 
systematic errors in the synthetic model spectra becomes pos- 
sible. 

A different meth od to tackl e this problem would be an 
approach like that of iDonati et aL C1994.) . who optimised di- 
rectly the frequency distribution of field strengths and direc- 
tions. Such an approach has the advantage that a formal best 
fit to the observed spectra is found which shows the smallest 
residuals of all possible combinations of database spectra, but 
a few potential traps should be kept in mind: (i) there is no 
unique relation between the B-\p diagram and the field topol- 
ogy (see Paper 1 for an example of ambiguous field configura- 
tions); (ii) there is no guarantee that the derived B-ij/ diagram 
corresponds to a field which fits in a globally consistent pic- 
ture if all rotational phases are regarded; (iii) and it cannot be 
guaranteed that a distribution of electric currents exists which 
produces the derived B-^ diagram. Remaining discrepancies 
between the observation s and integrated syn thetic spectra de- 
rived with the method of Donati et al. ( 1994) would be, never- 
theless, a good measure for the magnitude of systematic errors 
in the model spectra. 

Somewhat unfortunately, the two objects we have analysed 
so far with our Zeeman tomography code both have equal T^s 
(Paper 11 and this work). It would be desirable to examine the 
magnetic field geometries of MWDs for a broader range of ef- 
fective temperatures in order to search for potential effects of a 
secular field evolution as a function of cooling age. 

The outcome of the first three data releases of the Sloan 
Digital Sky Survey (SDSS) has nearly tripled the number 
of known MWDs. The multitude of newly found objects 
covers a broad range o f effective ternperatures and surface 
dipole field strengths ( Gansic ke et a l2Q02i: ISchmidtet alJ 
l2003t IVanlandingham et aL.200^ . Several new objects with 
B > 200 MG have been found, while the majority of objects 
is found in the low-field regime with B < 20 MG. The SDSS 
objects provide a vast hunting ground for further systematic 
studies of the field geometry of MWDs. The SDSS has also 
nearly doubled the number of helium-rich MWDs and an in- 
teresting option is the extension of our Zeeman tomography 
technique to the synthesis of their spectra. Calculations of the 
atomic parameters for He are available no w and fir st applica- 
tions to helium-rich MWDs are available djordan et aLi2001i: 
IWickramasinghe et alJ2002l) . 



Another promising avenue of research is the study of ac- 
creting MWDs in close mass-transferring binaries. In the fol- 
lowing papers of this series, we will investigate the magnetic 
field structures of the accreting MWDs in magnetic cataclysmic 
variables. 
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